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SUMMARY 

This study is devoted to a comparative analysis of three "Adaptive ZOOM" 
(ZOom Overlapping Multi-level) methods based on similar concepts of 

hierarchical multigrid local refinement : L.D.C. (Local Defect Correction), 

F.A.C. (Fast Adaptive Composite), and F.I.C. (Flux Interface Correction), 
which we proposed recently. These methods are tested on two examples of a 
bidimensional elliptic problem. We compare, for V-cycle procedures, the 

asymptotic evolution of the global error evaluated by discrete norms, the 

corresponding local errors, and the convergence rates of these algorithms. 

INTRODUCTION 


The need for local resolution in physical models occurs frequently in 

practice. Special local features of the operator coefficients, source terms, 

and boundary conditions can demand resolution in restricted regions of the 
domain that is much finer than the required global resolution. The multigrid 
methods with local mesh refinement provide one solution method to achieve 

efficient local resolution by solving problems on various locally nested 
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grids, and by using these grids as a basis for fast solution and correction on 

the global basic grid of the calculation domain. Different techniques have 

been proposed in the literature, such as the pioneering works [1, 2,3,4, 5]. 

Therefore, the concept of "Computational Adaptive Zoom" in the context of 

a "Graphical and Computational Architecture" has been introduced in the field 
of numerical simulation in order to take the best advantage of the new 

capabilities of high performance computer architectures [6]. It can be viewed 
as a generation made automatically (i.e. in an adaptive way) or not, of some 

multilevel hierarchical local nested zoom grids (ZG), overlapped all over the 

global basic grid (BG). These grids may move all over the entire computation 
domain Cl during the solution phase. This concept is supposed to allow both 
local refinement and global correction of the basic grid solution by a 

successive transfer of information between the connected grids (BG) and (ZG). 

So it is well adapted to a graphical vision of Zoom in terms of the creation 
of local graphical windows where it is needed in the problem (strong 

gradients, discontinuities, singularities,...), but in an active sense, i.e., 

the basic grid solution is modified and improved as the computing is 

performed. This has involved us in the creation of an original engineering 
software package called "AQUILON", still currently in development [6]. 

In addition, this strategy offers other interests. The goal is to combine 
the best features of both multigrid techniques and domain decomposition 

methods (in the case of overlapping grids) to provide an acceleration of the 
convergence rate and a good suitability for implementation on parallel 

computers, thus reducing the ellapse time. Moreover, another advantage is the 

possibilty to solve different differential problems on the grids (BG) and 

(ZG), which allows us to optimize both the physical and the numerical model. 

This can be particularly interesting for the approach of solving problems by 

"imbedding inside fictitious domains" associated with appropriate "control 

terms" for expressing the boundary conditions, as proposed in [6]. It is also 

possible to adopt different kinds of discretization on each grid. Thereby, the 

multigrid zoom methods share with the domain decomposition techniques the 

opportunity for obtaining precise solutions by combining solutions to problems 
posed on physical subdomains, or, more generally, by combining solutions to 
appropriately constructed continuous and discrete boundary value sub-problems. 

From the numerical point of view, the strategy adopted enables us to work 
only on structured and uniform meshes for each grid separately, on which a 

moderate number of degrees of freedom is required. On each grid, a "simple and 
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inexpensive" discretization is performed, leading to the same simple form of 
sparse pattern matrices (e.g. 2D block-tridiagonal). We aim at avoiding 

solving problems on unstructured or nonuniform composite meshes, which tends 

to introduce inaccuracies in the discretization, slowness in the solvers, and 

being surely more expensive in terms of implementation, data structures 
storage and CPU time. Our choice is expected to be relatively good in terms of 

duality quality/cost of computation for a lot of cases of moderate complexity. 

MULTIGRID ZOOM ALGORITHMS 

Different ZOOM algorithms will be examined and compared. We consider 
first the L.D.C. (Local Defect Correction) algorithm proposed by Hackbush [1]; 

we choose for the restriction operator a 2D bilinear interpolation one of type 

full weighting control volume". The second one belongs to the class of F.A.C. 
(Fast Adaptive Composite Grid) methods from McCormick [5], for which the 

analogy with the B.E.P.S. method [4] can be noticed. We use here the "delayed 

correction" version of F.A.C. Only the third one, the F.I.C. (Flux Interface 
Correction) algorithm that we proposed more recently [7], will be briefly 

described hereafter. 

All these Multigrid Zoom Algorithms are based on the same general 

principle : a successive transfer of information level by level, leading to 

the global correction of the initial discrete solution on each grid, and thus 
on the global basic grid (BG). The multilevel implementation is made in a 

recursive way as in the usual multigrid techniques (V-Cycles, W-Cycles, etc .) 

[1,3]. The resolution on each grid may be performed "exactly" or by using an 
inexact solve (e.g. a few iterations of a smoothing procedure). 

Notations and Definitions 

Consider the following second order non-linear elliptic boundary value 
problem defined on D a bounded, open domain in (R d , for d = 2 or 3 : 

(T) ( = div((p(u) ) + G(u) =f(x) x € Q (1) 

' well-posed boundary conditions on T = dQ symbolically called by (BC) 

The equation (1) L(u) s f is so expressed by splitting the nonlinear 
operator L(u) in the divergent part where cp(u) has the physical meaning of the 

flux density of the solution u= u(x) and the nonconservative one G = G(u). The 
relation between the solution u and the flux tp can take the general vector 
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form <p(u) = F(u) in many systems of conservation laws, but applications will 
concern an advection-diffusion equation or a Navier-Stokes problem. For the 

experiments here, (T) is a diffusion problem and we have (p(u) = -o.gradu . 

In order not to have too many formal requirements and restrictions, we 
assume explicitly only that this equation (1) has at least one isolated 

solution u* in the space L 2 (Q). All other assumptions are implicitly contained 

in the following considerations. 

The basic notations will be those classically used in the multigrid 

framework [1]. We denote by £ the current index of the grid level (0^ £ ^ £ ), 
{ = 0 is the level of the global basic grid (BG) which discretizes the entire 
calculation domain Q, and £ = £ ^*0 is the level of the most nested and finest 

zoom grid (ZG). Each grid of level £ can be characterized by : 



- the open domai n = | 0 j 

- the boundary = | | on which can be 

defined the unit outside normal vector 

- the c losure U 

- the mesh size h^ 


Each grid of level £ is 
associated to the nodes x e . 
successive grids of level £ and £+1 

be taken 
A„ = Q 


divided into a set 


We denote by r^ + j 
have Vl, 


and we 


successive mesh sizes will 
notations will also be used : 


l 


7+1 


as h t+1 = h ( / 
n and 7^ = H^ +1 


of control volumes V 

X 

the interface between two 
n n £+1 * 0. The 
2 P , p e IN*. The following 

n il. • 


The transfer operators between the grids t and £+1 

respectively, by R.jj , for the restriction operator and by 

prolongation operator. For all three algorithms, we have chosen P 


will be 

p£+l 

£ 


called, 
for the 


7+1 


F £,£+l ° U £ 


^£+ 1 ^ ^ 7+1 


£ 

n n 


as 


which is a monodimensional linear interpolation operator defined on the 
interface of the grids £ and £+1. Each value u £+1 at a node y e r £+l^ r £+l n ^ 
on the interface is obtained by a linear interpolation of the values u £ at the 
two neighbour nodes x et x’ belonging to (r^ ), and thus verifying 

Vi (y) = u £ (x) if y = x • 
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If we denote by u^ = the discretized equation of (1) on the grid of 
level £ , we can define the following discrete boundary value problems on 17^ : 


( P ) 
v o' 


L u = f : in Q 
ooo o 


on r 


: (BC) 




' L t u e - h ■ in n t 

- on n T : (BC) 

on \ (f { n f) : u { = P^_j 


£-1 


\ £*0 


( 2 ) 


k • __ 

We denote by u„ the discrete solution obtained on the grid 17. at the k-th 

1 k * k 

iteration of the zoom algorithm, and e^ = u^ - u^ the associated discrete 

error, where u^ is the natural restriction of the exact solution u* of problem 
(?) on n £ . 

For 0 < l < l , y(i) will represent the number of iterations of the zoom 

algorithm on the grid level £ in order to describe a whole cycle : if we have 
Y(£)=l (respectively y(£)=2), V0< £ <£ , then V-cycles (respectively W-cycles) 

will be described. We have y(l ) = 1, and y(0) is the total number of cycles 
performed from the basic grid (BG) in order to obtain the so-called 

convergence of the zoom algorithm. When £*=1 (i.e. for a two- grid algorithm), 

only V-cycles are of course carried out. The term "No Zoom" will be used for 
the resolution by "an exact solve" of problem (?) on the basic grid (BG) of 
mesh size h Q (1=0, k=0). The term "Zoom" will be used to indicate that some 
iterations of a multilevel zoom algorithm have been performed: t^Q, l<k<y(0). 


Description of The Multilevel F.I.C. Algorithm 

The main idea of the two-grid FIC method for levels i and £+1 is to give 

the opportunity to apply the local "flux residual" correction due to the whole 

patch level £+1, at each node x € ~K ^ on the grid level £. This is obtained 

through the expression of the local flux balance (i.e. integration) of eq.(l) 

over the volume n between the grid levels £, on one hand, and 

£+1 on the other hand. 

Because of the consistency of the conservative discretization of the 

fluxes by the finite volume method, which must be respected on each grid, the 

outside normal fluxes of (p(u) through an interface of two neighbour control 
volumes are opposite. By giving more importance to the local "flux residual", 
that leads to consider for the correction step on the grid level £, the local 

flux of the defect only at each node of a boundary zone 1^ £+ j defined as the 
"flux correction interface". We can choose for I^ + j , either the stripe A + = 

{u , x e d"Kji = n Q^}, or A' (see further Figure) if we want the 


boundary 9A to correspond to interfaces between control volumes on the grid 

level i : we will have V = V in the latter case. We define on the grid level 
t, r x u r, Vx e = l ll+l n , where V^= 9V^ n 9^ * 0, or 

respectively, r = 9V x n 3A^, , (mes (T^) = h^ ). 



We then propose the following restriction operator on the outside normal 
flux through the "interface boundary" £+1 = {u , x e 1^} : 


R £ + l ( W u) Vl )(x) 


R £+l : Y £,£+ 1 ° ^£+1 
1 


X £,£+l ° Q i 


mes(r ) 


Vi (u) -Vi dY Vx e ht+1 n q £ (3) 


We can then define, as in [7], the local "flux residual" correction at 
each node x e 1^ = I £g+1 n on the grid level £ by : 


r £ (<p)(x) = 


to(£,x)(u ) 

e(£,x) 


R £+l ( Vl (u) - n l+l> - ^Vl 


(x) 


(4) 


The control parameter e(£,x), which has the dimension of a length, has 
already been encountered in order to assign Neumann and Robin (or Fourier) 
boundary conditions in the context of "imbedding inside a fictitious domain" 
[6]. Its expression is given by : 


mes (V ) 

e(£,x) = (5) 

mes (T ) 


A complete calculation, still not published, gives a complex expression 

for co(£,x)(u), which is the following one in the case under subject of G = 0 : 
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jfc 

CO(£,x)(u) « 


r 

k * 

r 1 

(p(u).n dy 

r 

r J 

£+1 

f 

k 

p 1 

<p(u) . n dy 

r 

J T ' 

r 

k * 

(p(u).n dy 

r 

x ; 

f 

l+l - 

' 

(p(u).n dy 

r , J t 


We can then generate the successive iterates 
algorithm implemented in a recursive way : 


by the 


( 6 ) 


multilevel FIC 


Ini tialization : comp ute u ° 

u” is obtained by resolution of problem (7> o ) 

K 

Iterations : compute the successive iterates u Q 

for k = 1 to y(0 ) do FIC(O) 

Composite re - ac t ual i za t i on : providing uj (0) on (BG) by assigning 

for £ = £*-1 toOby step of -1 : uj ( 0 } (x ) = uJ+j } ( x) Vx e 

Pr o cedure FIC ( t ) 

If i=i* Then solve problem (7^*) Else 
b egin 

* 1 st step - resolution on the grid level t + 1 : 

- solve problem (7 > « + j) providing u f + j 

- for k = 1 to y( £+1 ) do FIC(£+1 ) 

* 2 nd step - correction on the grid level l : 

- solve problem (7^) with f ^ = f ^ + Xj r £(<p) 

where r^(cp) computed by equations (3)(4)(5)(6) 

and y t is the characteristic function of I p in II. 
e nd 1 i 

Remarks : 

1) - In any case, in order to avoid the explicit calculation of to(£,x)(u) 
by eq.(6), an economical solution is to use an approximate correction for FIC. 

In that version, called FIC(co), only the flux integrals on the interface 
will be evaluated by quadrature formulae (Simpson), and an average weighting 
factor co(£) will be determined by a semi-empirical way for each grid level. 
Besides, it can play the role of an average relaxation parameter for the 
iterative zoom algorithm when co(£)=co, . 
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2) - In terms of domain decomposition, the two-grid FIC for levels £ and 
£+1 can be regarded as a full overlapping iterative algorithm that splits the 

whole composite problem in two Dirichlet/ Neumann boundary value sub-problems: 

- the problem on the grid level £+1 with a Dirichlet boundary condition 

on the interface (2), 

- the problem on the level £ with a condition of relaxed transmission of 

the flux on the interface through (4), which demands the flux 

continuity at convergence. That condition can be considered as a Neumann 

boundary condition on y^ +1 by the technique of "fictitious domain" in [6], 

General Comments on the Three Algorithms 

i) - The two-grid FAC method for levels £ and £+1 can be regarded as an 

iterative procedure to solve "exactly" the discrete composite problem coming 

from an adequate discretization of problem (T) on the composite grid 

defined by the association of the grids and n^ + j. Therefore, the principle 

is to apply a multigrid algorithm between the grids U ^ and IZ on one hand, and 

between the grids and 17^ ^ on the other [5,4], There is therefore a 

correction phase on both the grid levels £ and £+1 with respect to the 

discretization on the composite grid. In that sense, FAC can be viewed as an 

"exact" solver for the composite problem. Because the composite grid stencils 

agree with the coarse and fine grid stencils, respectively, outside and inside 
the refinement region, and because the correction equations are solved 

exactly, the composite grid residual is nonzero only at the interface. 

ii) - Due to the attention needed for the nonuniform discretization of the 

problem on the interface zone of the composite grid, FAC method can prove to 
be a little difficult to implement in a more than two grids version. 

iii) - On the contrary, LDC and FIC methods, which are easier implementing 

in the multilevel case, are only approximate solvers : they don’t use a 

composite grid and neglect the fine grid residual correction. The former 
consists in the local correction of the solution defect inside A £ as the 

latter involves a local flux residual correction through the interface Y £ £ + j- 

iv) - Both FAC and FIC methods provide corrections by balancing fluxes 
computed from both coarse and fine grids across the interface. They take the 

best advantage of a conservative discretization of the equations, for example, 
by a finite volume technique. 


NUMERICAL APPLICATIONS 


In that context, we propose to compare three types of multigrid zoom 

algorithms on two examples of a linear elliptic problem (T) presenting, 

respectively, a discontinuity of the operator coefficients for (Tl) [8], and a 
singularity of the exact solution for (T2) [1] : 


(?) 


L(u) = -d i v( a( x ) .gradu) + a(x) u = f(x) 

- a, a >0 e L°°(Q) et f g L 2 (Q) 

well -posed boundary conditions on T = dQ 


in Q =]0,1[ X ]0,1[ 
symbol i cally ca 1 led by 


(1 ’) 
(BC) 


These problems were already tested successfully on the FIC method in [7]. 

Problem (Tl) is heterogeneous and defined by M), a^O, a^lOO inside a disk of 
radius 0.1 and o^l outside (Fig. la). A solution computed on a very fine basic 
mesh (512 2 ) will be used as the reference solution u . Problem (T2) is defined 

by M), a^O, osl (Fig. lb); the exact solution is u* = ln(r) with r= V x 2 + y 2 . 


Numerical Implementation and Procedures 


The discretization on each grid, independant of the geometry of the 

problem, is made in a conservative way by a finite volume method on a uniform 

Cartesian mesh. The classical five-point scheme is used providing a second 

order precision. The resolution of the linear systems, which are 
block-tridiagonal and symmetric positive definite, is performed by a fast and 
efficient solver : a preconditionned conjugate gradient (PCG) method (CG-SSOR) 

vectorized by a Red and Black numbering of unknowns. The results for two grids 

are obtained by an "exact" solve on each grid. The results for multilevel LDC 

or FIC (£ > 1) are given for an "inexact" solve on each grid (including 

Fig.5b), i.e., a fixed number itcg of iterations of PCG on each grid with : 

itcg = 2 for h Q = 1/8 itcg = 4 for h Q = 1/16 //eg =8 for h Q = 1/32 

The results are analyzed with different norms (L°°, L 2 , L-energy norm) of 

the discrete error evaluated on the global basic grid (BG, £=0). We study the 

asymptotic evolution of the relative error norms I le”l I / I ki Q l I (No Zoom) and 

^Y(°) _ ||ej (0) || / | lu j | (after y(0) Zoom iterations) with e^ = u^ - u Q , which 
allows us to estimate an asymptotic average rate x : 
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for l = 1, as function of h ] or p (for a fixed h Q ) 

o ' 1/m 


X = 


kV 0 , (p^) 


, wi th m = max 1 p e IN 


| P G IN* | 


Here m=3 and y(0)=2, see Tab.l, Tab.2, and Fig.2a, Fig.3, and Fig.4. 

a f i > 

\ 1/m 


* for i > 1, as function of i (for a fixed h Q and p=l) 


X = 


“n 




, wi th m = max 


| l e IN* | 


Here m=3 and y(0)=10, see Tab. 3, and Fig. 2b. 


The convergence rate of LDC, FAC and FIC have been also compared (Tab.4): 

* for FAC : we study the variations of the Euclidean norm of the 

v 

composite residual | r 0 ( u )| 2 for k = 1 to y(0) (Fig.5a), and a convergence rate 


p is calculated by a geometric mean : 


P = 


' i ^ <0> (“)i 2 ' 

. I rj (»>!, . 


i/(y(o>-i) 


Ic 1c 

* for LDC or FIC: we study the variations of quantities 5 q = \ |u Q 
for k = 1 to y(0) (Fig.5b), and a convergence rate p is then estimated by : 


u o _ 1 iL 2 


P = 


gY(o) ^ l/(y(o)-l) 

u o 


si 


Comparative Numerical Results 

1) - By comparing a no-zoom method and a ZOOM one, we notice that the 

error globally decreases ; between two increments of p or l , it is divided by 

an elevated average t-factor of between 1.5 and 3.5 (Tab.l, Tab.2, Tab.3). For 

problem (T2), the decrease is monotonic and there seems to be good analogy 

* * 

between the variation of the error as a function of p (for l =1) or of l (for 
p=l), (see Fig. 2a and Fig. 2b). The influence of the position and dimensions of 
the local grids (ZG) becomes negligible as h Q decreases [7], Due to the choice 
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of discretizing on a Cartesian mesh independently of the geometry of the 

heterogeneity, the error for problem (T\) does not decrease monotonically as 
already noticed in [7,8]. 

2) - In many cases, the error obtained with zoom is less than computed 

without zoom on a single basic grid of mesh size h < h * In particular, Fig.4 

k oi 

shows that the local discrete error |e^| calculated point by point on the 
diagonal of the domain (?2), by a two-grid FIC method (h Q = 1/1 6, h^hg/2, k=2) 

is globally better than the error obtained with No Zoom at the corresponding 

nodes of BG (1=0, h Q = 1/32). The former results are more accurate inside the 
refinement region and get closer to the latter case far from the interface. Such 

remarks can also be made for the discrete error norms in the other Tables or 
Figures. However, the error is not reduced beyond a threshold value consistent 
with the order of precision of the discretization schemes on the different 

grids (cf, the multigrid defect correction method using Richardson 
extrapolation [1]). 

3) - The two-grid FAC and FIC methods yield error results of the same 

order of magnitude for both problems. These results are far better than for 

LDC for problem (?>1), where the flux conservation plays an important role. On 
the contrary, LDC yields as good results as the others for problem (P2), and 

sometimes better. However, as LDC does not deal with the interface fluxes, but 

only works on the solution inside the open refinement region, it can become 
inefficient (t = 1) if the refinement region does not contain enough coarse 

nodes on which the local defect correction is performed (Tab.l, Tab.2, Tab.3). 

4) - The results with the version FlC(to) for 0,1 < 0) < 0.5 are nearly 

similar to those obtained with co = co(£,x)(u) calculated by (6) (Fig. 3). That 

could justify the interest of the approximate version FIC(co), and particularly 

as a preconditioner of the discrete problem, as suggested in [4], 

5) - Because of its exact character, the FAC method yields the far best 

convergence rate, a mean value of 0.16, nearly independant of both h Q and h ] 

(Fig.5a and Tab.4). We obtain a mean convergence rate of 0.42 for FIC(co=0.2), 
just a little better than LDC with 0.50 . These convergence rates remain not 

very sensitive to the variations of h Q and £ (Fig. 5b and Tab.4). However, 

those of FIC have a noticeable tendency to become better as the number of grid 
levels (or i) increases (see Tab.4). 
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CONCLUSION 


Despite its non-exact character, FIC provides as good results as FAC, 
concerning the analysis of discrete errors for both the two tested problems. 
In particular, FAC and FIC proved to be better than LDC for problems where the 
flux conservation locally plays a main role. 

FAC yields very good convergence rates (p-0.16), better than LDC (p-0.50) 
or FIC (pO.42), but its multilevel implementation remains more difficult. 
However, the use of FIC as a preconditioning technique of the discrete problem 
is likely to be very interesting, especially for the approximate version 

FIC(O)) where the factor co becomes a relaxation parameter. We are currently 

testing such a procedure for Navier-Stokes problems. 
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Tab. 1. 


Problem CPI) - Two Grid Zoom ?(0) = 2 - Discrete L norm of the error 


h o 

NO ZOOM 

h 

v* 

ZOOM x^O.375 et x a =0. 625 

ZOOM x^O.25 et x z «0. 75 

LDC 

FAC 

FIC 

u * 0.55 

LDC 

FAC 

FIC 

w = 0.40 

1/8 

0.434E-1 

p=l 

0. 434E-1 

0. 493E-1 

0.210E-1 

0. 434E-1 

0. 732E-2 

0.639E-2 

p=2 

0. 434E-1 

0. 790E-1 

0. 122E-1 

0.434E-1 

0. 133E-1 

0. 106E-1 

p=3 

0.434E-1 

0. 330E-1 

0. 589E-2 

0.434E-1 

0.289E-2 

0. 251E-2 

T 

1.00 

1.10 

1.95 

1.00 

2.47 

2.59 

1/16 

0.689E-2 

P=i 

0. 689E-2 

0. 167E-1 

0. 119E-1 

0. 739E-2 

0. 986E-2 

0. 106E-1 

p=2 

0. 689E-2 

0. 374E-2 

0. 394E-2 

0. 297E-2 

0. 123E-2 

0. 161E-2 

p=3 

0.689E-2 

0.255E-2 

0. 339E-2 

0.214E-2 

0. 673E-3 

0.592E-3 

T 

1.00 

1.39 

1.27 

1.48 

2.17 

2.27 

1/32 

0. 982E-2 

p=l 

0. 136E-1 

0. 342E-2 

0. 538E-2 

0.244E-2 

0. 160E-2 

0.201E-2 

p=2 

0. 140E-1 

0.295E-2 

0. 412E-2 

0. 174E-2 

0. 428E-3 

0. 392E-3 

p=3 

0. 134E-1 

0. 326E-2 

0. 331 E-2 

0. 153E-2 

0. 233E-3 

0. 262E-3 

T 

0.90 

1.44 

1.44 

1.86 

3.48 

3.35 

1/64 

0. 192E-2 


1/128 

0.648E-3 

1/256 

0. 193E-3 


Tab. 2. Problem (P2) - Two Grid Zoom r(O) = 2 - Discrete L-Energy norm of the error 


h 

0 

NO ZOOM 

h 

h 

1 2 P 

ZOOM x^O et x 2 =0.25 

ZOOM x t =0 et x 2 *0. 5 

LDC 

FAC 

FIC 

w * 0. 20 

LDC 

FAC 

o° 

= ! 
3 

1/8 

0. 342E-1 

p=l 

0.342E-1 

0. 152E-1 

0. 168E-1 

0. 120E-1 

0. 138E-1 

0. 139E-1 

p=2 

0. 342E-1 

0.817E-2 

0. 110E-1 

0. 441E-2 

0. 511E-2 

0. 553E-2 

p=3 

0. 342E-1 

0.622E-2 

0.990E-2 

0.221E-2 

0. 250E-2 

0. 346E-2 

T 

1.00 

1.77 

1.51 

2.49 

2.39 

2.15 

1/16 

0. 206E-1 

P =1 

0.723E-2 

0. 829E-2 

0. 837E-2 

0. 704E-2 

0. 818E-2 

0.819E-2 

p=2 

0.270E-2 

0. 308E-2 

0. 337E-2 

0. 232E-2 

0. 282E-2 

0. 285E-2 

P-3 

0. 140E-2 

0. 152E-2 

0. 215E-2 

0. 627E-3 

0. 975E-3 

0. 107E-2 

T 

2.45 

2,38 

2. 12 

3.20 

2.76 

2.68 

1/32 

0.133E-1 

p=l 

0. 454E-2 

0. 527E-2 

0.527E-2 

0. 453E-2 

0. 526E-2 

0. 526E-2 

p=2 

0. 150E-2 

0. 182E-2 

0. 184E-2 

0. 149E-2 

0. 180E-2 

0. 180E-2 

P-3 

0. 407E-3 

0.628E-3 

0.699E-3 

0.387E-3 

0.570E-3 

0. 577E-3 

T 

3.20 

2.77 

2.67 

3.25 

2.86 

2.85 

1/64 

0.888E-2 


1/128 

0.607E-2 

1/256 

0. 420E-2 


1/512 


0.294E-2 




Tab. 3 


Tab. 4. 


Problem (T2) - Multilevel Zoom LDC/FIC 

a'(O) = 10, j = h^/2 0 ^ l s £ -1 , = 0 and x = 0.5 - 

2 

Discrete L norm of the error 


h 

HO ZOOM 

number 

of 

ZOOM 

ZOOM FIC 

0 

• 

l = 0 

grids 

LDC 

to - 0. 2 

to * 0.35 



l m = 1 

0. 806E-2 

0. 869E-2 

0. 787E-2 

1/8 

0.209E-1 

/ - 2 

0. 440E-2 

0.625E-2 

0. 465E-2 

- 3 

0. 368E-2 

0. 606E-2 

0. 428E-2 



T 

1.78 

1.51 

1.70 



- 1 

0. 399E-2 

0. 404E-2 

0. 394E-2 

1/16 

0. 105E-1 

l* - 2 

0. 157E-2 

0. 199E-2 

0. 166E-2 

e* = 3 

0. 102E-2 

0.165E-2 

0. 115E-2 



T 

2.18 

1.85 

2.09 



** = 1 

0. 186E-2 

0 . 202E-2 

0. 249E-2 

1/32 


i* * 2 

0.969E-3 

0. 768E-3 

0. 809E-3 

U. ocVt-Z 

** = 3 

0. 71 IE- 3 

0. 416E-3 

0.514E-3 



T 

1.95 

2.33 

2.18 

1/64 

0.265E-2 





1/128 

0. 132E-2 





1/256 

0.662E-3 





1/512 

0. 331E-3 






Problem (T 2) - Multilevel Zoom LDC/FIC - Two Grid Zoom FAC 
y(0) = 10, h^ +J = h^/2 p 0 s 1 s ! -1 , x ( = O and x z = 0.5 - 
Convergence rate p 


h o 

number 

of 

grids 

ZOOM 

LDC 

ZOOM FIC 
w = 0.2 

ZOOM 

FAC 





p = 1 

p = 1 

p=l 

p=2 

p=3 


• 

i 

* 

1 

0. 53 

0.46 

0. 14 

0. 15 

0. 15 

1/8 

c 

K 

2 

0. 45 

0.42 




t* 


3 

0.43 

0.34 





~ LJ ' • 

l 


1 

0.56 

0.44 

0. 14 

0. 16 

0. 16 


l * 

* 

2 

0.49 

0. 42 




1/ lo 

l 9 


3 

0. 48 

0.39 





m 

1 

I 

0.50 

0.47 

0. 15 

0. 17 

0. 18 


,/32 







No Zoom 
Zoom 1/h • B 


Zoom 1/h » 16 
0 

-Eh Zoom 1/h * 32 
0 










